Optimal geometry for efficient loading of an optical dipole trap 
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One important factor which determines efficiency of loading cold atoms into an optical dipole trap 
from a magneto-optical trap is the distance between the trap centers. By studying this efficiency for 
various optical trap depths (2-110 mK) we find that for optimum dipole trap loading, longitudinal 
displacements up to 15 mm are necessary. An explanation for this observation is presented and 
compared with other work and a simple analytical formula is derived for the optimum distance 
between the trap centers. 
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I. INTRODUCTION 

Optical dipole trapping is an established tool for storing, manipulating and studying cold atomic gases P, 0, fl- 
optical dipole traps (ODTs) helped to achieve Bose-Einstein condensate of several elements [1, H, 0, 0, Q as well as to 
make ultra-stable frequency standards ODTs are also widely used for studies of various phenomena at quantum 
degeneracy, such as the Mott insulator - superfluid liquid transition [lo| . 

The simplest realization of an optical dipole trap (ODT) is a single, tightly focused Gaussian laser beam, tuned far 
below the atomic resonance frequency. The trapping potential is proportional to the light intensity which, for a single 
beam, results in a highly elongated trap shape with optical potential minimum at the trapping beam focus. 

The ODT needs to be loaded by a precooled gas, and a standard way of loading atoms into the ODT is to transfer 
them from a magneto-optical trap (MOT). Kuppens et al. [IT] discuss ways of maximizing this transfer by optimizing 
various parameters: MOT Hght intensities and detunings, magnetic field gradient, the ODT depth, and alignment of 
the traps. In particular, they found that controlling the geometry of the overlap of the MOT and ODT beams allows 
for substantial improvement of the transfer efficiency. They found that the loading rate is optimum with a longitudinal 
displacement between the centers of the ODT and the MOT, and that this optimum displacement increases with the 
ODT depth. Optimal displacements between the ODT and MOT centers reported in Ref. were about one-half 
of a MOT diameter (the depths of the studied ODTs were 1-6 mK). These observations were rather qualitative and 
not accompanied by systematic quantitative analysis. 

Since the loading of the trap is a crucial step in any experimental work involving ODT, the purpose of our work was 
to determine quantitatively the optimum displacement between the ODT and MOT centers for maximum transfer of 
atoms in a wide range of the ODT depths. The trap depths were changed by varying the ODT laser power from 50 to 
660 mW and detuning from 0.43 to 1.72 nm. The experimental setup and procedure are described in the next Section, 



while the obtained results are presented in Section UlI 
results with the models described by Kuppens et al. 



and in Section HVl we compare quantitatively our experimental 
iH and O'Hara et al. 



II. EXPERIMENT 



We use a standard MOT setup with three pairs of orthogonal, retroreflected beams for trapping *^Rb atoms. The 
trapping light, red detuned from the 5^Si/2 F = 2 — > 5^P3/2 f = 3 transition, has a total six-beam power of 15 mW. 
The repumping light (0.5 mW, tuned to the 5^Si/2i^ = 1 ^ 5^P3/2i^' ~ 2 transition) is added to one pair of the 
trapping beams. The radius of the beams is 8 mm (1/e^ intensity). The resulting MOT radius is 0.66 mm (1/e^ 
density). After collection, the atoms are further cooled by increasing the detuning of the MOT beams to 100 MHz 
and decreasing the repumper intensity (see Fig. [T]). At the end of the MOT cooling stage the MOT temperature 
reaches Tmot = 17 /iK, as determined from free expansion of the released atom cloud. 
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FIG. 1: The time sequence of loading the ODT from the MOT. 

The dipole trap is formed by a focused beam of a Ti:Sapphire laser. The ODT beam is first expanded, then focused 
by a 25-cm focal length achromatic lens to a waist of 12 /im. The corresponding Rayleigh length is = 0.54 mm. 
We apply the ODT beam with powers, Pqdt, of 50, 100, 200, 400, and 660 mW, and red detunings from the D2 line, 
AAoDT, of 0.43, 0.86, and 1.72 nm. The resulting trap depths determined by the ratio Podt/AAqdt were J/o/fce 
= 2.1, 4.2, 8.4, 17, 33, 67, and 111 mK. For the trap-depth estimation, a possible effect of the Rb Dl line has been 
neglected as it does not exceed a few percent. 

The transfer of the Rb atoms from the MOT to the ODT proceeded with the timing shown in Fig. [U it is similar 
to the sequence used by Kuppens et al. [ll| . 

For each ODT laser power and detuning, the ODT-focus position was gradually varied relative to the MOT center. 
The described loading sequence (Fig. [1]) was repeated for each position of the traps and the number of atoms loaded 
into the ODT was monitored by integrating the fluorescence image captured by the CCD camera. 

III. RESULTS 

Figure [2] shows the results of the loading measurements performed for the trap depth of J/o/fce = 33 mK. The 
fluorescence images of the dipole trap are presented in Fig. [2] (a) for several different positions z of the ODT relative 
to the MOT center. Superimposed are the contours of the original MOT (1/e^ density contour) and the dipole trap 
(the equipotential line at U — —2.bkBTuoT)- The diffused clouds seen below the trap contours are the atoms which 
were not captured into the ODT, falling in the direction of gravity after MOT has been switched off. Figure [2] (b) 
presents the total fluorescence from the ODT, proportional to the number of trapped atoms, versus the focus position. 
It is evident that when the centers of MOT and ODT coincide (ODT focus position z = 0), the number of captured 
atoms is over 5 times smaller than in the case of displaced ODT: z = ±10 mm. The observed asymmetry between 
positive and negative z positions is attributed to geometrical imperfections. We believe that for perfect alignment of 
the laser beams and magnetic field, the I{z) dependences in Fig.[2](b) should be symmetric. The difference in heights 
of the two curves shown in Fig. [2] (b) is most likely caused by some drift of the MOT conditions. 

We have conducted the described measurement for a range of ODT trap depths, J/o/fce, between 2.1 and 111 mK. 
For the trap depth of 2.1 mK, we obtain a single, broad loading maximum at z = 0. The width of this maximum 
(FWHM) is 9 mm. For increasing ODT potential depth, the maximum broadens and eventually, at 8.4 mK, splits 
into two separate maxima such as shown in Fig. [2] (b). 

With our fixed beam waist, the optimal separation of MOT and ODT depends only on the ratio Podt/AAqdt, 
that is, it depends only on the ODT depth. We have verified this for trap depths of 8.4, 17, and 33 mK realized by 
different ODT laser powers and detunings. We have verified that optimum displacements z corresponding to a given 
-Pqdt/AAodt ratio were constant within the experimental error of a few percent for various combinations of Pqdt 
and AAodt which preserved the Podt/AAqdt ratio. Figure [3] presents the dependence of the optimum displacement 
on the trap depths. It is worth noting that in the case of the deepest optical traps realized with our maximum ODT 
beam power, this optimum separation reaches 15 mm. This displacement is over 10 times larger than the diameter 
of the MOT, whereas the optimum displacements reported in Ref. [l3| were only of the order of one-half of a MOT 
diameter. 
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FIG. 2: Loading of the 33-mK-deep dipole trap for different positions of the ODP laser focus with respect to the MOT center, 
(a) Fluorescence images of atoms in the ODT, vertically stretched for better visualization of the ODT potential. The elliptic 
contours mark the MOT positions (0.66 mm radius) before switching on the ODT. The elongated contours represent the ODT 
equipotential surfaces, U = — 2.5fcBTMOT. (b) Total fluorescence from the ODT, proportional to the number of atoms loaded 
into the ODT, versus the focus position. Solid line: Podt = 200 mW, AA = 0.43 nm; dashed line: Podt = 400 mW, AA = 
0.86 nm. 



IV. DISCUSSION 

The dipole potential depends on the local light intensity /(r) and detuning from the atomic transition 5X as 

2^2c3 r 



where T denotes transition Hnewidth [ll| . The intensity of a focused Gaussian beam of power P is described by 

2P 



I{r,z) 



■ cxp —2 



where (r, z) are the cyHndrical coordinates, and: 



(1) 



(2) 



(3) 



is the 1/e^ beam radius. The characteristic dimension in the radial direction is the beam waist wq, while in the axial 
direction it is the Rayleigh length zr = ttwq/X. 

The resulting trapping potential is visualized in Fig. [H The potential is harmonic near the focus, with ellipsoidal 
equipotential surfaces. However, the more distant equipotential surfaces are qualitatively different from the inner 
ones, and acquire characteristic peanut shape. From equations ^ and ^ one can deduce that the equipotential 
surfaces change their character when U becomes smaller than U' = Uo/e, where Uq is the trap depth and e is the 
base of the natural logarithm. The equipotential surfaces U = const with \U\ < \U'\ have two locations Zmax where 
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FIG. 3: Optimum displacement of the optical dipole trap relative to the magneto-optical trap for maximum transfer of atoms 
versus the ODT depth. The black dots represent our experimental data. The solid line is calculated from Eq. ([5| with the 
scaling parameter a = 2.5, as described in the text. The dashed line is the prediction based on the model of O'Hara et al. with 
no free parameters [12]. See also Sec. IIVI 
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FIG. 4: The cross sections of the equipotential surfaces of the ODT formed by a focused Gaussian beam, z/zr is the distance 
along the beam in units of the Rayleigh length, r/wo is the radial coordinate in units of the beam waist. Uo is the trap depth. 
Dashed lines mark the positions of izmax for U = Q.lUo- 



their circumference is maximized: 



Zinax = ±^R\/--F7 ^ 1- (4) 

V e (7 

We will now compare our experimental data with the two models of the ODT loading described in the literature 

Kuppens et al. pjj] assume that the initial loading rate is proportional to the MOT atomic density niviOT, the 
mean velocity of atoms w, the effective surface area of ODT, A, and the trapping probability, Ptrap- The first three 
parameters contribute to the flux of atoms into the ODT volume. When studying geometrical effects, it is important 
to consider the effective surface area A. As the Rayleigh length is much larger that the size of the MOT, the radii of 
the ODT equipotential surfaces change insignificantly over the MOT region (see also Fig. [SKa) below). The effective 
surface area across which atoms are loaded into the ODT is then proportional to the local radius of the ODT. The 
crucial question in this reasoning is: which of the equipotential surfaces should represent the "surface of the ODT". 
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Kuppens et al. assume the surface U ~ kBTuoT, where the ODT potential U is comparable to the kinetic energy of 
atoms in the MOT. 

In order to compare quantitatively our data with the model described above, we write the equation describing the 
ODT surface as U — afceTMOT and look for a scaling factor a, of the order of unity, which gives the best fit to our 
data. Equation ^ now takes the form 



^ = ±W--|^-1. (5) 
V e afcsJMOT 

Using the least-squares method, we find that Eq. |[5|) best describes our experiment for a = 2.5. The result of this 
fitting is shown by the solid line in Fig. [3l It turns out that despite its simplicity, the model by Kuppens et al. gives 
a good qualitative description of our data, but needs the scaling factor a to be determined from the experiment. 

Another related work is that of O'Hara et al. [l3| who developed a microscopic model of ODT loading based on the 
Fokker-Planck equation. Their model is capable of describing the spatial and temporal changes of the atom density 
during ODT loading. For comparison with our results, the relevant equation is the expression for the initial number 
of atoms iV(0) contained in the trap Eqs. 11, 12) 



2 fl 

7V(0) = (TT^/^iJno) Inw) exp(-gw), (6) 

2 Jo 

where w is the beam radius (the width of the ODT potential in the radial direction), R and uq are the radius and 
the density of the MOT cloud, respectively, and q = C/o/fceTMOT- The calculations of O'Hara et al. were performed 
for the case when MOT and ODT centers coincide, but remain valid also if the centers are displaced axially. Then 
Uq in Eq. ^ is treated not as a global trap depth Uq = mm{U{r, z)}, but as a local trap depth for a particular axial 

r,z 

position z: 

C/o(.) = minM.,z)} = ^,^^ (7) 

(compare Eqs. HI) and On the other hand, the local beam radius w(z) increases from its initial value wo as one 
moves away from the focus, according to Eq. ([3]). Consequently, as the axial coordinate z increases, the trap becomes 
shallower and more energetic atoms remain untrapped. On the other hand, the trap also becomes broader and atoms 
are captured from a larger volume. The balance of these two competing mechanisms results in a maximum loading 
for a certain displacement between the MOT and ODT centers. 

The initial number of trapped atoms has been calculated according to Eq. ^ for different displacements z, treating 
the parameters w and q as z-dependent. The resulting function is plotted in Fig.[5]for the trap depth [/q/^b — 33 mK. 
Note the maximum at z = ±13 mm, close to our experimental result: z — ±10 mm (see FiglS]). Similarly, we calculated 
the maxima of the N{z) function for different trap depths and the result is shown by the dashed fine in Fig.O It agrees 
qualitatively with the experimental data, but for deeper traps the predicted optimum displacement is overestimated. 
Nevertheless, the agreement is remarkable for a model that does not require any free parameters. It should be stressed 
that the model described by O'Hara et al. was developed for lithium atoms in CO2 laser traps, where the trap-induced 
light shift of the atomic transition frequency is small fl^. This is not the case of our experiment, hence this may be 
the main reason of the discrepancy between our data and the model of O'Hara et al. 

Under conditions of our experiment, the frequency shifts for the MOT trapping and repumping transitions may be 
taken with good accuracy as 2U{r,z)/h, where U{r,z) is the ODT potential which takes into account solely the D2 
line. As it was already mentioned above, under our experimental conditions, the excited state makes the most 

important contribution to the ODT potential and light shifts. Other states, and ^£'3/2,5/2) contribute nearly 

two orders of magnitude less and can be neglected at the current level of experimental precision. The light shift also 
modifies the trapping transition but, given the large detuning (100 MHz) of the MOT beams in the loading phase, 
its effect on the trapping-beam detuning is negligible. On the other hand, since the repumper beam is resonant, the 
repumping efficiency at 2; = Zmax,?" = is reduced by about a factor of four (for all trap depths) which allows to 
accumulate atoms in the dark F = 1 ground state. In Fig. [6l we compare the spatial ranges in which the light shift 
generates the repumper detuning with the spatial MOT extension calculated for the trap depth of 33 mK. As can be 
seen in Figure Eljb), this reduction affects only the central part of the MOT which, consequently, results in a "dark 
spot" effect which is known to increase the MOT density [iJl and enhance the trap loading process [l3|- Figure [6ljc) 
shows that the light-shift-induced detuning of the repumper beam, and consequently the dark spot effect, decrease 
for increasing z. For z = 13, which would be the optimum displacement predicted from the model of O'Hara et al, 
the repumper detuning amounts only 2.5 MHz and the dark spot effect is much weaker. Thus, the fight shift and 
the related dark-spot effect decrease the optimum MOT-ODT displacement relative to the predictions of Ref. [l^ . 
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FIG. 5: The initial number of trapped atoms, A'^, as a function of the MOT-ODT displacement, z, predicted by the model of 
O'Hara et al. [l^ for our experiment (trap depth = 33 mK). 
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FIG. 6: The geometrical superposition of the trap potentials and the influence of the ODT beam on the MOT. (a) The 
geometry of the MOT (represented by the l/e'^ density outline) and the 33-mK-deep ODT (represented by the equipotential 
surface U = — 2.5feBTMOT - see text). The displacement of the traps is optimized for maximum atom transfer, (b) The 
detuning of the MOT-repumper caused by the presence of the ODT beam, plotted along the dashed line in (a). The gray 
outline represents the atomic density distribution of the MOT. (c) Same as (b), but along the ODT beam axis (r = 0). Note 
the vertical scale has been expanded for better visualization. 



as illustrated in Figure [31 We believe this is the main source of discrepancy between our results and the model of 
O'Hara et al. which refers to the case of negligible light-shifts. 

Another factor of potential importance for the ODT dynamics is the photon scattering rate. With a trap relatively 
close to resonance, such as ours, the spontaneous photon scattering leads to atom heating and, eventually, to atom 
loss. We have estimated the heating rates under conditions of our experiment. For the deepest trap of 111 mK 
the scattering rate is 430 kHz and the heating during the 50 ms holding time is estimated to be about 2 mK. This 
estimation was done for the trap center, whereas, as we demonstrated, the atoms are most efficiently trapped far off 
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the center (15mm), where the Hght intensity is significantly lower and the local scattering rate is only 530 Hz. Thus, 
we believe that it is safe to neglect this mechanism in the case of our experiment. 



We have performed a systematic study of the efficiency of the ODT loading from a MOT as a function of geometrical 
arrangement of the ODT and MOT in a wide range of ODT depths (2 to 110 mK). We have observed that for the 
optimum loading the trap centers need to be displaced in the direction of the ODT beam. While similar conclusion 
was formulated earlier in Ref. [llj, we were able to verify it in a more than 20-times wider range. We confirmed that 
the optimum displacement depends only on the trap depth, i.e., on the ratio between the trap Hght power and its 
detuning. Analyzing the optimum displacements, we were able to quantitatively compare our results with the existing 
models of Kuppens et al. [ll| and O'Hara et al. [I4I. Our data are well described by the former model after fitting of 
a scaling factor and agree qualitatively with the latter model without any free parameters. 

Based on the present work, we propose the following semi empirical formula for the optimal trap displacement: 
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